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Abstract. We examine under which conditions one may 
apply, to steady state keplerian accretion discs, the 0- 
viscosity prescription which has been derived from rotat- 
ing shear flow experiments (v = pflR 2 , where f2 is the 
angular velocity at radius R and f3 is a constant of order 
10~ 5 (Richard & Zahn 1999)). Using a vertically aver- 
aged model, we show that this law may be suitable for all 
three families of known systems: in young stellar objects, 
evolved binary stars and Active Galactic Nuclei discs (ex- 
cept in their outer gas pressure dominated regions where 
turbulence becomes hypersonic). According to the stan- 
dard criterion for viscous stability, /3-discs are always sta- 
ble throughout. Using realistic opacities and equation of 
state, we demonstrate that these discs are thermally un- 
stable in the temperature domain where hydrogen recom- 
bines, when they are optically thick, and this could lead 
to limit cycle behavior. Radiation pressure dominated re- 
gions are thermally stable, in contrast with a-discs. This 
results in a fully stable solution for the innermost parts of 
AGN discs. 

Key words: accretion, accretion discs — instabilities — 
hydrodynamics — galaxies: active — galaxies: nuclei 



1. Introduction 

What causes the transport of matter and angular momen- 
tum in accretion discs? This question is the subject of a 
longstanding debate, but most of those working in this 
field tend now to agree that it has been answered, since 
even a very weak magnetic field leads to a linear insta- 
bility (Chandrasekhar 1961) which, in its fully developed 
regime, is able to provide the necessary stresses (Balbus & 
Hawley 1991; Hawley & Balbus 1991). However, it is un- 
likely that this mechanism operates in cold neutral discs 
(e.g. Gammie & Menou 1998). Besides, the issue which 
has not been settled yet is whether this magnetic instabil- 
ity is the only instability responsible for the "anomalous 
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viscosity" , or whether a purely hydrodynamic instability, 
generated by shear in differentially rotating discs, may 
not be of comparable or even higher efficiency. Since a 
non-magnetic keplerian disc is linearly stable, this hydro- 
dynamic instability needs a finite amplitude perturbation 
to be triggered and laboratory experiments indicate that 
it appears only above Reynolds numbers of the order of 
10 4 . Such conditions are still out of reach of the most am- 
bitious numerical simulations; therefore, these cannot be 
used yet to prove or to disprove the occurrence of shear 
instability. 

It is worth noting that the widely used a-prescription 
for turbulent viscosity (Shakura & Sunyaev 1973) does 
not invoke any particular mechanism at the origin of tur- 
bulence: it is a simple parameterization which is tailored 
to yield turbulent velocities that remain subsonic (or sub- 
alfvenic) for a < 1. It is often asserted that the main 
success of the a-prescription is the prediction of recurrent 
outbursts of dwarf novae (e.g. Cannizzo 1993). However, 
if such eruptive phenomena are indeed of thermal origin, 
as is commonly thought, it is likely that many ad-hoc vis- 
cosities, not only the a-prescription, could work as well. 

For these multiple reasons it makes sense to investigate 
the basic properties of accretion discs built with the alter- 
nate /3-prescription, which is observed in laboratory rotat- 
ing shear flows, as explained below. This is our intention 
in this paper. We first present in Sect. 2 some properties 
of /3- viscosity, and examine in Sect. 3 to which regimes it 
may be applied, allowing for subsonic turbulent velocities 
in geometrically thin, radiatively cooled accretion discs. 
In Sect. 4 we study the viscous and thermal stabilities of 
/3-discs, using standard criteria, and draw our conclusions 
in the last section. The Appendix contains the vertically 
averaged equations for /3-discs and the derivation of the 
criterion of thermal instability. 

2. Comparing the a and /3-viscosities 

The laboratory experiment which is the most pertinent 
to studies of hydrodynamical instabilities in differentially 
rotating flows is the Couette- Taylor experiment: a fluid is 
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sheared between two cylinders rotating at different speeds. 
Only a few experiments have been run in the case where 
the angular momentum increases outwards, as in a kep- 
lcrian disc, but from the torque measurements that are 
available, it appears that the turbulent viscosity scales as 



v t oc R 3 



dR 



(1) 



where R is the distance from the rotation axis and £1 is 
the angular velocity (Richard & Zahn 1999). In a keplerian 
disc, this relation is equivalent to 



up = priR 2 



(2) 



which we shall call the /^-prescription. 

This prescription was originally proposed by Lyndcn- 
Bcll & Pringle (1974) and some aspects of it have been dis- 
cussed in later papers (De Freitas Pacheco & Steiner 1976; 
Thompson et al. 1977; Lin & Papaloizou 1980; Williams 
1980; Hubeny 1990). It has been recently revived by Bier- 
mann and Duschl (1998; see also Duschl, Biermann & 
Strittmatter 2000) on the grounds that the only relevant 
scales in a keplerian disc are the angular velocity and the 
radius, since they contain all the information about the ro- 
tation and the curvature of the flow. These authors make 
the sensible choice of equating the parameter (3 with the 
inverse of the critical Reynolds number IZe, which they as- 
sume to be the same as in plane-parallel shear flows (i.e. 
Tie 10 3 and therefore (3 » 10~ 3 ). 

However, in the following we shall use a smaller value 
for this /3-parameter, which is derived from the Couette- 
Taylor experiment. The only experimental data available, 
in the case where angular momentum increases outwards, 
are those obtained by Wendt (1933) and Taylor (1936), 
in which the inner cylinder is at rest; from these one de- 
rives (3 ~ 10~ 5 . This value is compatible with the critical 
gradient Reynolds number which characterizes these ex- 
periments (see Richard & Zahn 1999): 



Tie* 



R? 



dn 



dR 



6 x 10 5 



(3) 



where v is the molecular viscosity of the fluid. Clearly, 
new experiments are required to measure this parameter 
in regimes which resemble more the keplerian flow; in the 
meanwhile, until new results are obtained, we shall adopt 
the value quoted above. 

One important property of this /3-viscosity is that it 
depends only on the radius in a keplerian disc, and does 
not involve local physical conditions. On the contrary, the 
standard a-prescription (Shakura & Sunyaev 1973) de- 
pends on the local values of the pressure scale height H 
and sound velocity c s (taken in a broad sense - it is the 
Alfvcn velocity if magnetic pressure dominates): 



v a = ac s H. 



(4) 



To emphasize the difference between the present prescrip- 
tion and the standard law, we use the relation c s = ClH 



which expresses the hydrostatic equilibrium in a keplerian 
vertically averaged disc, and rewrite Eq. ([!]) as 



H 

v n = aflR I — 
R 



" ;X v U 



(5) 



We see that the ratio v a jvp depends on the aspect ratio 
of the a-disc. In general, the quantity (H/R) a is weakly 
sensitive to both R and a, but it depends strongly on the 
central mass and the accretion rate, and is typically in the 
range 0.001 - 0.1 (e.g. Hure 1998). 

The a-prescription has been designed such that the 
Mach number is smaller than unity for a < 1. In con- 
trast, subsonic turbulence is not guaranteed with the (3- 
prescription. This restriction on subsonic turbulent veloc- 
ities is dictated by the fact that supersonic turbulence 
would be highly dissipative and therefore difficult to sus- 
tain. Moreover, the turbulent viscosity used here has been 
measured in a liquid, and its application to compressible 
fluids can only be justified in the subsonic range. There- 
fore one has to check whether the /3-prescription predicts 
subsonic turbulence in a given disc. By making the reason- 
able assumption that the vorticity ft/^t of the turbulent 
eddies is of order Q, (v t and it denoting respectively the 
typical velocity and length scales of turbulence), and set- 
ting vp = Vtlt, the Mach number in a /3-disc is given by 



c s c s V R 



(6) 



where (H/R)p is the aspect ratio of the /3-disc. We note 
that A4 t < 1 provided H/R > yfj3, which is not a too 
severe requirement with our low value of /3 (the situation 
would be more problematic with /3 = 10~ 3 , conflicting 
with the keplerian assumption). 

3. Application range of the /3-viscosity 

Let us now examine under which conditions the (3- 
prescription allows for subsonic turbulence in geometri- 
cally thin accretion discs in active galactic nuclei (AGN), 
evolved binary systems (EBS) and young stellar objects 
(YSO). For this first exploration, we will be content with 
the usual one-dimensional steady state approach to com- 
pute the sound speed and the aspect ratio, the vertical 
structure being reduced to vertical averages. The reader 
can find the relevant equation set in Appendix A. 

In a stationary keplerian disc characterized by a steady 
accretion rate M, the surface density E is determined by 
the conservation of mass and momentum: 



M = 3 7w£, 



(7) 



ignoring a correcting factor which is negligible sufficiently 
far from the center (see Frank et al. 1992). Applying the 
/3-prescription (0), we find 
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where f ec id — 4.5 M M s 1 is the Eddington factor 
(around 0.1 typically), M = M 8 x 10 s M Q = M M , 
M = M M Q = mi 6 x 10 16 g/s = M_ 8 x 10~ 8 M /yr, 
i? = as i? s = Rio x 10 10 cm = r AU AU (i? s = 2GM/c 2 
being the Schwarzschild radius of the black hole), and 
/3_g = (3 x 10 5 . These dimensionless variables are well 
suited to describe the three families of discs quoted above. 
Note that Eq. (JsJ) is quite general and does not depend on 
the vertical structure. 

Thermal equilibrium between radiative cooling and 
viscous heating provides another equation. Assuming that 
the disc is optically thick, we have (Frank, King & Raine 
1992): 

16crT 4 3Q 2 M 



3k r E 8tt ^ 

where a is the Stefan constant and kr the Rosseland mean 
opacity. 

3.1. Regimes compatible with subsonic turbulence 

In gas pressure supported discs, the hydrostatic equation 
yields directly the mid-plane temperature T: 



(10) 



p A*m H 

with the usual notations for the gas pressure P g , the den- 
sity p, the Boltzmann constant k, the molecular weight p 
and the mass of the hydrogen atom mn- It then follows 
from Eqs. (§), d), d) and © that the turbulent Mach 
number in those discs is 



M t 



3.1 /rj s 

0.3 /3% 8 
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V*Ml'*Ml'"T-W. (11) 



Note the relatively strong dependence of the Mach number 
on the viscosity parameter and its weak dependence on the 
opacity and specially the radius. 

We conclude that the /^-prescription is certainly appli- 
cable to gas pressure dominated accretion discs around 
forming stars, white dwarfs, neutrons stars and stellar 
black holes, but not to such discs around supermassive 
black holes, where this prescription would imply super- 
sonic turbulence. 

The innermost part of accretion discs may be domi- 
nated by radiation pressure if they reach a sufficiently high 
temperature; this occurs for instance in standard discs sur- 
rounding black holes or compact objects. In such a case, 
the hydrostatic equation is written 



p 3c E 



(12) 
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Fig. 1. Disc thickness versus radius for a /3-disc with (3 — 
10~ 5 (solid line) and for a standard disc with a = 0.01 
(dashed line), for = 1 M Q and M = 10~ 8 M /yr. 

and therefore, according to (^|), H depends only on the 
accretion rate, and not on the viscosity prescription: 



H 



3k r 



M. 



16-7TC 

Thus the Mach number is given by 



(13) 



3 x 10- J /3 
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dd ■ 
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I Mq X 



(14) 



meaning that turbulence is expected to be subsonic in the 
inner part of radiation pressure dominated discs. 

3.2. Thin keplerian discs assumption 

In this paper, we restrict ourselves to keplerian discs, 
which are necessarily geometrically thin, i.e. H/R <C 1. 
Eq. (^) shows that the turbulent Mach number increases 
as the aspect ratio decreases; therefore the disc cannot be 
thinner than y/~j3 x R. One easily finds from Eq. (Illl) that 



1.5 x 10 



-1/2 -7/16 • 1/4 1/16 

5 A* 1 m r 10 



0.08 4^ £3 V 1/9 mo 7/16 r^ (15) 



within the gas pressure dominated regions. It is interesting 
again to make the comparison with standard discs. Using 
power law solutions of a-discs (Hure 1998) and Eq. (|l5|), 
we obtain 



f^0.9/3„ 5 



l/8 a l/10 



(16) 



indicating that /3-discs globally should have about the 
same shape as standard a-discs. This is confirmed in 
Fig. U, where we have plotted the disc thickness computed 
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from vertically averaged equations with realistic opacities 
and equations of state, for the two viscosity laws and for 
M = 1 M and M = 1CT 8 M Q /yr which could repre- 
sent both a YSO disc and a EBS disc. In the innermost 
part of AGN, where the radiation pressure dominates, the 
disc thickness is independent of the viscosity prescription, 
as we have already pointed out, and therefore Hp = H a . 
We thus conclude that, despite the flaring, /3-discs are ex- 
pected to be thin and keplerian, provided that the accre- 
tion rate does not reach too large values. Note that a-discs 
are subject to the same constraint. 

4. Stability of /3-discs 

We shall now examine two kinds of instability that an 
accretion disc can undergo: viscous instability which can 
lead to disc fragmentation and thermal instability. 



4-1- Viscous stability 

A general condition for viscous stability is (Frank et al. 
1992) 



d\uM 
<91n£ 



> 0. 



(17) 



In a standard a-disc, the equilibrium relation M(S) takes 
the famous "S" -shape with the instability occuring on the 
intermediate branch of negative slope (Cannizzo 1993). 
On the contrary, according to Eqs. (||) and (fjj) /3-discs 
are viscously stable everywhere, since Ineq. (|l7j ) is always 
fulfiled. We illustrate this important property in Fig. ||, 
with two typical M(E)-relations obtained through vertical 
structure computations for a /3-disc and a a-disc (Hure, 
2000). 

4-2. Thermal stability 

Thermal instability occurs when the energy balance be- 
tween radiative losses and viscous heating can no longer 
be maintained after perturbation of the equilibrium state. 
The condition for stability is written (Pringle 1981): 



dlnQ+ 
d\nT 



d\nQ- 
d\nT 



< 



(18) 



where Q + is the heating rate (right-hand-side of Eq. (g)) 
and Q~ is the cooling rate (with Q + = Q~ at equilib- 
rium). The above criterion is quite general and can be 
applied to non standard viscosity prescriptions, provided 
that the vertical dynamical time scale t z ~ H/c s is smaller 
than the thermal time scale tth ~ PH/Q + . The require- 
ment of a subsonic turbulence discussed in the previous 
section automatically ensures that t z < t t h since we have 
the relation 



tth 

t, 



1 



H 



1 

A4T 



(5-disc 
a-disc 



log t (in g/cm ) 

Fig. 2. Relation between the accretion rate and the sur- 
face density, from vertical structure computations, at 
R = 10 10 cm from the center and with M=l M©, us- 
ing the a-viscosity prescription (with a = 0.1) and the 
/3-prescription (/3 — 10 -5 ). 

It is easy to see that Q + cx f2 3 i? 2 S with the /3- 
prescription and so the thermal stability of a /3-disc is 
only governed by properties of the cooling function. The 
situation is different in an a-disc where the heating rate 
Q + cx f2 3 _ff 2 E is implicitly a function of T through H. 
To determine the slope of the cooling function Q~(T), we 
have used vertically averaged equations (see Appendix A) 
and computed the term (91nQ~/<91nT) analytically (see 
Appendix B). 

The resulting expression Eq. ( |F3 . 7| ) may be used to 
analyze the thermal stability for very simple ideal cool- 
ing functions: electron scattering and free-free processes, 
as considered already by Piran (1978). We have distin- 
guished 4 cases: optically thick, optically thin, gas pres- 
sure dominated and radiation pressure dominated. This 



(19) 



can easily be handled by suitable settings in Eq. (B.7) (see 
Eqs. (B.8) and Eq. (B.ll)). The results are summarized in 
Tab.([lj). The remarkable point is that /3-discs are always 
thermally stable for these ideal cooling mechanisms, which 
is in complete agreement with the conclusions drawn by 
Piran (1978). Note that the optically thin gas pressure 
dominated regime with free-free emission is marginally 
stable whereas it is unstable with the standard viscosity 
law, as noticed by Pringle, Rees & Pacholczyk (1973). 

The situation is quite different when using realistic 
Rosseland and Planck opacities as well as an accurate 
equation of state. A typical result is displayed in Fig. |[ 
where we have considered a disc surrounding a one solar 
mass central object and M = 10 -8 M /yr. We find that 
the disc is thermally unstable when the midplane temper- 
ature is in the range 5000 — 9000 K (zone 2 on the figure), 
due to the steep rise in the opacity. Like in a standard disc, 
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Table 1. Quantity (dln<3~ /91nT)s,n and inferred thermal stability for ideal opacity sources (disregarding self- 
gravitation). Values and comments in parenthesis concern the a-prescription. Various cases are considered: /3 = 1 (/3 
is the ratio of the gas pressure to the total pressure); (3 — \ and (3 — (radiative pressure only); Thompson scattering 
only (k = (Tx); free-free opacity only (n = /% oc pT~ 7 / 2 ); both absorption processes (with x = — ^ — < 1). 
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Fig. 3. Right panel: Temperature radial profile (thin line) 
and thermally instable regions (bold line) in a "realistic" 
/3-disc for M — 1 M© and M = 10" 8 M /yr. Left panel: 
the quantity (91nQ — /SlnT)!;^ versus T computed from 
Eq. (§[7|). 



Fig. 4. Same legend as for Fig. || but for M — 10 -12 
M0/yr typical of a quiescent EBS disc. The optically 
thick/thin limit is reached at no — 0.2 where T ~ 6600 
K. The optical thickness r goes through a minimum for 
T ~ 4500 K where r ~ 4 x 10" 4 . 



this unstable regime corresponds to the partial ionization 
of hydrogen. A specially interesting case is displayed in 
Fig. [| where the accretion rate is M — 10~ 12 M Q /yr , 
typical for an EBS disc in quiescence. There is a (multi- 
valued) optically thin solution with T < 6600 K (zone 
2bis in the figure). We have checked that this regime is 
also compatible with both H/R -C 1 and M.t < 1- We see 
that the hottest branch of this solution, which has already 
been investigated by Williams (1980), is thermally stable. 

Similar computations have been carried out for many 
input pairs (M, M) corresponding to AGN, EBS and YSO 
discs and the overall conclusions are the following: 

— the hot optically thick solution which extends down to 
the central object (zone 1 in Figs. 3 and 4) is always 



thermally stable (this includes the radiative pressure 
dominated solution); 

— the region where hydrogen recombines (zone 2) is al- 
ways unstable if the disc is optically thick; 

— the hottest branch of the optically thin solution with 
7000 K > T > 4000 K and dT/dR < (zone 2bis), 
when it exists, is always stable, even if hydrogen re- 
combination occurs; 

— the cold branches of the optically thin solution with 
T < 4000 K and dT/dR > are always unstable; 

— the outermost parts (zone 3) are globally stable. Thin 
unstable zones can be present. 
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5. Conclusion 

In this paper, we have examined some properties of thin 
accretion discs built with the /3-prescription for the turbu- 
lent viscosity, which has been derived from the transport 
of angular momentum observed in differentially rotating 
laboratory flows (Richard & Zahn 1999). This prescription 
may be applied whenever it predicts a turbulent velocity 
that is subsonic, and we have verified that, with the value 
we have taken for (3 (i.e. 10 -5 ), this condition is fulfhed in 
most discs, around stellar objects or massive black holes, 
except in the outer part of AGN discs where gas pressure 
dominates. Only radiatively cooled discs have been sur- 
veyed; it would be interesting to investigate ADAF-type 
solutions with this prescription. Let us recall that our con- 
clusions have been obtained here with vertically averaged 
models, and it is plausible that a more realistic 2D-model 
would extend somewhat the domain of applicability of the 
/3-prescription (Hure & Richard 2000). 

Although the /3-viscosity and the a-viscosity may have 
comparable magnitudes, they have very different effects 
both on the structure and on the stability of steady ke- 
plerian discs. An important property of /3-discs is that 
they are viscously stable and do not tend to fragment 
into concentric rings, contrary to a-discs. Also, unlike a- 
discs, they are thermally stable for ideal cooling processes, 
such as Thomson scattering and free-free absorption, as 
was shown by Piran (1978). However, with more realistic 
opacities and equation of state, models of /3-discs contain 
a thermally unstable region corresponding to the recombi- 
nation of hydrogen, much like standard a-discs. It remains 
to be seen with time-dependent models whether this in- 
stability can lead, as we suspect, to a limit cycle behavior. 

Acknowledgements. The authors wish to thank the referee for 
his competent remarks and suggestions. 
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Appendix A: Vertically averaged equations for a 
steady state keplerian /3-disc 

In a steady state keplerian accretion disc where the viscos- 
ity is given by Eq. @) , the mass accretion rate M is linked 
to the surface density S by the relation 



M = 3 7r/3fLR 2 £, 



(A.I) 



where f2 is the keplerian angular velocity and R the radius. 
The surface density is £ = 2pH, with p being the mass 
density and H the total pressure scale height, which is 
determined from the equation of hydrostatic equilibrium 



(A.2) 



where ( represents the effect of self-gravity. The total pres- 
sure P is the sum of the gas pressure P g and of the radi- 
ation pressure P r : 



P = P K + P r 



We draw P g from the perfect gas equation 
pkT 



prriH 



(A.3) 



(AT) 



where T is the midplane temperature and p the mean 
mass per particle in units of the proton mass. The disc 
temperature in the equatorial plane is fixed by the balance 
Q + = Q~ between viscous heating (Frank, King & Raine, 
1992) 



+ _ 3Q 2 M 
and radiative losses (Hubeny 1990) 

8crT 4 



Q~ = 



(A.5) 



(A.6) 
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where kr, and up are respectively the Rosseland and 
Planck mean opacities (in cm 2 .g _1 ). The general expres- 
sion for the radiative pressure is then 



8 
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Note that in the optically thick limit the radiation pressure tends to its LTE value (labeled with an asterisk) 



4oT 4 
p — = p* 

r 3c 



(A.8) 



The above equation set can in fact be reduced to a system of two non linear algebraic equations with £ as the 

(A.9) 



unknown. From Eqs.( A.1)-(A.7), the first equation is 
„ 2 7r 6 /3 6 G 4 M 4 E 8 16 crT 4 



P* - ttGE" 

256 M 6 P 9ck r kp5(E) 
where the function <S(£) is 
2£ 4 



0. 



y; 2 



2 \/3kr 3k r k p 



The second equation is obtained from Eqs.( |A.l| ), ( |A.5|) and 



(A.10) 



(A.ll) 



We can follow the method described in Hure (1998) to solve the system of Eqs.( A.9 ) and ( A.ll ). Note that the 
solution is valid in all cases: optically thick or thin, pressure dominated by gas or radiation, self-gravitating or not. 

A.l. The optically thick limit 



If we set <S(E) = |S in Eqs.(A.E) and (ATI), and omit the last term of Eq.(A.9) we recover the optically thick limit 



(r ^> 1), and there is a single equation to solve: 

M 12/T + AiAf 2/7 + A = 0. 

The coefficients Ao and Ai are respectively 

6/7 



Ao = 16( 3Vgk r ) — and Al = -I^^ W?- 



2/7 



(A.12) 



(A.13) 



Note that the solution of Eq.(A.12) must satisfy t>1 



A. 2. The optically thin limit 



In the optically thin limit (t <C 1), S(T.) = 4/3krkp and the last term of Eq.(AT)) is identical to P* . There is also a 
single equation to solve, namely 



6/5' 



M 



3kp 2 /67r 5 /? b G 4 M 



T 3/5 



5fl6^4 M 4\ 2 /5 



CTKp 



(AAA) 



and its solution must be such that t«1. 



Appendix B: Criterion for the thermal instability 

Setting r = kr£ + 8/3 KpS + 4/\/3, and differentiating Eq.( [A.6|) while keeping £ constant, we obtain 



dlnQ- 
<91nT 



80 P 
3 Tp 



3 rp 



9 In p 
<91nT 



where tr = krE and rp = KpS are the total Rosseland and Planck optical thicknesses respectively, and 
' d In Kr \ / 9 in Kp \ / <9 In Kp \ / 9 In ftp 



or 



<91np 



6r = 



<91nT 



a P = 



9 hip 



6 P = 



91nT 



(B.l) 
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Note that this expression does not depend on the actual viscosity prescription. For Thompson scattering, we have 
simply ap = 6r = ap = bp = 0, and for free- free opacity, ap = ap = 1 and bp = bp = — Z. For a combination of both 
sources, we have ap = ap = x and bp = bp = — |x where x = Ke/(Kg + ctt)- 

To compute the derivative (<91np/<91nT) s n , we use the two different expressions of the total pressure at the 
midplanc. First, by differentiation of Eq.(jAT2|), keeping E and fl constant, we have 



<91nP 
SlnT 



1 



E.S2 



dhip 
1 + C \dlnT 



Note that (<91np/i91nT)n e is always negative. Second, from Eq.(A.3) which is written 
and setting xf = (dlnPg/dlnT) , xf- = (<91nP g /<91np) T and — P s /P, we have 



(B.2) 



(B.3) 



<91nP 
d\nT 



e.si 



/?x| + 1-/3 or 



TP 



73 J 



91np 
91nT 



+/3 X |+ 1-/3 



TR 



TR + 



&r+ 1-/3 



Combining now Eqs.(B.2) and (B.4), we find 



d\np 
dlnT 



i + C 



+ /3x| + (l-/9)«R- 



TR 



TR 



%/3 



^ + (1-/3) 



TR 



TR 



-bn+ (1-/3 



d\nQ- 
dlnT 



<91nT 



(B.4) 



E,Q 



(B.5) 



E.o 



Finally, from Eqs.(B.l) and (JB . 5|) and after rearrangement 



d\np 
dlnT 



4(1 - /?) - (1 - $) 



m&R - §£ 



+ ^ + (1-/3) 



TR 



(1-/3) 



8 ap 

3 Tp 



T ^+/3x|+(l-/3") 



TR + 7^ 
Tr 



or- 



Tr + ^/J 



(B.6) 



and so Eq.(B.l) becomes 
01nQ-\ 

E,0 



TP 



-bp + 



TR. 



i + C 



+ /3 X |+ 1-/3 a R 



Tr 



'4r-r R 6 R +|^ 



+ /3x|+ 1-/3 or 



TR 



T R 



+ 13-1 



V3, 



tr ir - 



8 ap 

3 Tp 



V3, 



(B.7) 



It is then easy to derive from Eq.(B.7) the following asymptotic expressions: 
in the optically thick limit (r ~ Tp and Tp > tr 3> 1) 



dlnQ- 
ainT 



3 x 8 + 4(l-3) 

= bR - 4 - QR 3 v g 7 1 > 

E,0 PXt + T+C 



— in the radiative pressure dominated limit (0 ~ 0) 
'91nQ- 



<91nT 



6 R -4-4a R (l + C), 



— in the gas pressure dominated limit (/3 ~ 1) 



aing- 



= Or — 4 — Or , 

2 4- r 

E,n z ^ s 



ainT 

in the optically thin limit (t ~ and Tp -C 1) 



'd\nQ- 



M + 4(1 + P)h 



= —bp — 4 + ap- 
E,n (1 - /3)ap + /3Xt + TTc 



(B.6 



(B.9) 



(B.10) 



(B.ll) 



